Modulated pair condensate of p-orbital ultracold fermions 
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We show that an interesting of pairing occurs for spin-imbalanced Fermi gases under a specific 
experimental condition — the spin up and spin down Fermi levels lying within the px and s orbital 
bands of an optical lattice, respectively. The pairs condense at a finite momentum equal to the 
sum of the two Fermi momenta of spin up and spin down fermions and form a p-orbital pair 
condensate. This 2kF momentum dependence has been seen before in the spin- and charge- density 
waves, but it differs from the usual p-wave superfluids such as ^He, where the orbital symmetry 
refers to the relative motion within each pair. Our conclusion is based on the density matrix 
renormalization group analysis for the one-dimensional (ID) system and mean-field theory for the 
quasi-lD system. The phase diagram of the quasi-lD system is calculated, showing that the p- 
orbital pair condensate occurs in a wide range of fiUings. In the strongly attractive limit, the system 
realizes an unconventional BEG beyond Feynman's no-node theorem. The possible experimental 
signatures of this phase in molecule projection experiment are discussed. 

PACS numbers: 03.75.Ss, 71.10.Fd, 37.10.Jk, 05.30.Fk 



I. INTRODUCTION 

Pairing with mismatched Fermi surfaces has long fas- 
cinated researchers in the fields of heavy fermion and or- 
ganic superconductors, color superconductivity in quark 
matter [l|, and, most recently, ultracold Fermi gases 
with spin imbalance In a- classic two-component 

model for superconductivity, the mismatch arises from 
the spin polarization of fermions in the same energy 
band. Its effect was predicted to produce intriguing, un- 
conventional superfluids such as the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) deformed Fermi surface HQ, 

and breached pair phases [l^, [HI . The limiting case of 
large spin imbalance was also studied to explore the for- 
mation of Fermi polarons [l2j . In parallel, the behavior 
of particles in the higher orbital bands of optical lattices, 
due to large filling factors, thermal excitations or strong 
interactions, is widely studied for novel orbital orderings 
of both bosons (Tsl - fisj and fermions [l^, with repul- 
sive interactions. Recently, interband pairing of unpo- 
larized fermions was shown theoretically to give rise to 
Cooper pair density waves [l^ . 

In this article, we report a fermion pairing phase re- 
sulting from the interplay of Fermi surface mismatch and 
p-orbital band physics. In such a phase, the pair con- 
densate wave function is spatially modulated and has a 
p-wave symmetry. This phase arises in an attractive two- 
component Fermi gas on anisotropic optical lattices un- 
der a previously unexplored condition of spin imbalance. 
Namely the majority (t) spin and the minority (J,) spin 
occupy up to Fermi levels lying in the Px and s bands, 
respectively. We show that pairings take place near the 
respective Fermi surfaces of the spin '|" fermions in px 
band and J, fermions in s band. This induces a modu- 



lated p-orbital pair condensate that differs from the usual 
p-wave superfluids such as '^He. The state requires only 
an on-site isotropic contact interaction and the pair is 
a spin singlet, while the ^He p-wave superconductivity 
has to involve anisotropic interaction and spin triplet. 
The modulation wave vector of the order parameter is 
Q w kp-f + kpi, where kp^f, kp^ are Fermi momenta for 
spin t and \, species, respectively. This 2kp momentum 
dependence is an unprecedented signature in superflu- 
ids other than the spin- and charge- density waves. In 
the strongly attractive limit, tightly bounded pairs con- 
dense at flnitc momentum Q, which realizes an unconven- 
tional Bose-Einstcin condensate beyond Feynman's no- 
node theorem [13i ,15i . JQ.. i2^1. 



II. MODEL 

The system under consideration is at zero tempera- 
ture and consists of two-component fermions in a three- 
dimensional (3D) cubic optical lattice with lattice con- 
stant a, described by the Hamiltonian 

rf'xV4(x)[-|^V^ + y(x) - M.]^,(x) 

+g I d3x^,|(x)Vj(x)V'i(x)V't(x). (1) 

Here "00- (x) is the fermionic fleld operator at x with spin 
a =t, i, V{x.) is the lattice potential, is the chemical 
potential for spin a fermions, and 5 < is the contact at- 
traction which can be tuned by the Feshbach resonance. 
In particular, we consider the case where the lattice po- 
tential in the x (parallel) direction is much weaker than 
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the other two (transverse) directions, so the system be- 
haves quasi-one-dimensionally. 

We expand i'ai'x) = J2nr <^n(x-r)c„r, where 0„(x-r) 
is the nth band Wannier function at lattice site r with 
Cnr the annihilation operator in Wannier basis. As a re- 
sult, we obtain the usual attractive Hubbard model with 
nearest-neighbor hopping between ith site with orbital 
band a and jth site with orbital band /3 



2m 



y(x) 



(2) 

■,)0,,(x-rj). 



and on-site attraction between orbitals 

Uafs.y,j =9 y (i^x(/)*(x-ri)0^(x-ri)0^(x-i 

(3) 

The lowest two energy bands are the s and Px band 
(the Py and Pz band are much higher in energy because 
of tighter confinement in the transverse directions). For 
brevity the Px band is simply called p band in the fol- 
lowing. By filling fermions with spin t to the p band and 
spin J, to the s band, the Hamiltonian becomes 



c.) - fis y n 



- J2 itsSlSr^' + t^PjPr" + h.C.) - /.p ^ < 
{r,r") r 

-Hw6^< + C/,p^nX- (4) 

r r 

Here, (r, r') and (r,r"} denote the nearest neighboring 
lattice sites in parallel and transverse directions, tl 
and tl are the hopping amplitudes along the parallel 
direction for the s- and p-band fermions respectively, 
while = tp = t'^ are the hopping amplitudes in 
transverse directions. (Pr) is the annihilation oper- 
ator at lattice site r for s-band J, (p-band t) fermions. 

= S^SrjTiP — P^Pr are the number operators, and 
fJ-s,lJ-p are the corresponding chemical potentials. Usp 
is the attractive on-site interaction between s- and p- 
band fermions and can be tuned by changing the scat- 
tering length using Feshbach resonance, wt, is related to 
the band gap. In the tight binding region we assume 
Wf, ^ \Usp\, and consequently the s-band fully filled spin 
'I" fermions are dynamically inert and not included in the 



III. DMRG CALCULATION FOR ID CASE 

First we consider the pairing problem in the simplest 
case of ID {t^ = 0), which is schematically shown in 
Fig. [ija). The two relevant Fermi momenta are kp^ 
(for s-band J, fermions) and kp^ (for p-band t fermions). 
From a weak coupling point of view, to pair fermions of 
opposite spin near their respective Fermi surfaces, the 
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FIG. 1. (Color online) (a) A schematic illustration showing 
the pairing between s- and p-band fermions. The s band is 
also fully occupied with f fermions (not shown), (b) The spa- 
tial variation of the pairing correlation C{x) for Ns = 49, 
Np — 15 according to DMRG. The blue scatters are the 
DMRG result and the solid line is the fitting using function 
a cos{qx + b)/x'^ + c. The inset in (b) shows the s- and p-wave 
Wannier functions in momentum space, which are elongated 
in the transverse direction (in real space they are compressed 
in transverse direction). The s-wave Wannier function has 
even parity while the p-wave Wannier function has odd par- 
ity. 



Cooper pairs have to carry finite center-of-mass momen- 
tum (CMM) due to Fermi surface mismatch. Further- 
more, in order for all Cooper pairs to have roughly the 
same CMM, the only choice is to pair fermions of opposite 
chirality. Note that the dispersion of p band is inverted 
with respect to the s band, so pairing occurs between 
fermions with momenta of the same sign but opposite 
group velocities. These elementary considerations show 
that the CMM of the pair should be approximately the 
sum of two Fermi momenta. 
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Fi 



(5) 



This result differs from that of the usual one-dimensional 
spin imbalanced fermions within the same band, where 
the FFLO pair momentum is the difference, Q ~ \kp-f — 
kpi I , as found in a two- leg-ladder system [2l| . 

Mean-field theory and weak coupling consideration can 
provide only a qualitative picture for ID problems. To 
unambiguously identify the nature of the ground state, 
we use density matrix renormalization group (DMRG) 
to compute the pair correlation function. In the numer- 
ical calculations, we used parameters tl = 1 as the unit 
of energy, tf, = 8, /i^ = 1.7, fj,p — ujb = ^ H, hr which 
the ratio between t^ and tp is chosen according to typ- 
ical tight-binding bandwidth ratio. Usp is tunable with 
Feshbach resonance and in the following calculation we 
will focus on Usp ~ —9 p^ . The truncation error is con- 
trolled in the order of 10~^ or less. Equation ^ predicts 
Q Ri kp^ + kpi = 0.4357r/a. Figure [Ijb) shows the pair- 
ing correlation function in real space Cij = {SjpJ PjSj) 
as a function of a:: = \ for a chain of = 60 sites with 
open boundary condition, where the indices i and j are 
real space positions. Since the system only has algebraic 
order, C{x) decays with x according to a power law. On 
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top of this, however, there is also an obvious oscillation. 
A curve fit with formula C{x) = acos{qx + b)/x^ + c, 
shown in Fig. [Hb), yields a period q ~ 0.4387r/a, 
which is in good agreement with the wave number given 
by Eq. ([5]) before. The Fourier transform of the pair cor- 
relation function 



^9 - TV ^ - 



(6) 



is peaked at q = 0.4267r/a (to be plotted in Sec. IVip . 
These features of the pair correlation function are the 
signature of the existence of the 2kF CMM pairing in 
our system [H, H] • 



IV. MEANFIELD ANALYSIS FOR QUASI ID 
CASE 

Now we move on to the quasi- ID system where a weak 
transverse hopping t-^ <^ is added. We carry out a 
mean-field analysis of Hamiltonian Hgp by introducing 
the s-p pairing order parameter 



Ar — Usp{SrPr) J 



(7) 



where (...) means the ground-state expectation value. 
Two different trial ground states are investigated, the ex- 
ponential wave Aj. = Ae**^ "", which is analogous to the 
Fulde-Ferrcll phase and the cosine wave Ar = A cos Q • r, 
which is analogous to the Larkin-Ovchinnikov phase. Q 
and A are determined self-consistently by minimization 
of ground-state free energy (i?sp)- Transverse hopping 
introduces a small Fermi surface curvature and spoils 
the perfect nesting condition as in the pure ID prob- 
lem above. However, the curvature is small for weak t-^. 
Thus, we expect Q pointing almost along the parallel di- 
rection, Q = (3(1,0,0), in order to maximize the phase 
space of pairing. 

The mean-field Hamiltonian for the exponential wave 
can be diagonalized in momentum space by standard pro- 
cedure. We get the ground state energy 



(Hsp) 



k,7=± 



sp 



with the self-consistent gap equation for A 
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Here, k is lattice momentum, N'^ is the total number of 
sites, 6 is a step function, and Aj^'''' = ^[^k ^ ^Q-k ^ 



4A2 -I- (^^ -I- ^Q_i(.)^] eigenenergy of the Bogoli- 

ubov quasiparticles. As evident from these formulas, the 
pairing occurs between an s-band fermion of momentum 
k and a p-band fermion of momentum Q — k with dis- 
persion — —2tl cos k^a — 2t^ cos kyU — 2t^ cos k^a — /is 



and = 2tp cos kxa ^ 2t'^ cos kya^2t-^ cos k^a — fip + ujh, 
respectively. 

The cosine wave is spatially inhomogeneous. A 
full mean-field analysis requires solving the Bogoliubov- 
de Gennes equation to determine the gap profile self- 
consistently. Here we are interested only in computing 
the free energy for the ansatz Ar = A cos Q • r to com- 
pare with the exponential wave case. Thus, it is sufficient 
to numerically diagonalize the full Hamiltonian Eq. ([4]) 
for a finite size lattice. We introduce a vector of dimen- 
sion 2N 



(SI 



.s 



k^ , ky, O ^ 

where A:" = 2Tm/Na is the discrete momentum 
in the x direction. The components of a obey 
anticommutation relation a^^^]} — (5miin2i 

where TOi,m2 labels the corresponding operator com- 
ponent of a. The Hamiltonian takes the compact 



form Hsp = J2 



Sq^-q)N^A'^ /2Usp. Since Hkyk^ is real and symmetric, 
it can be diagonalized by an orthogonal transformation 
akyk, = 'Dkyk^Pkyk, to ylcM 27V eigenvalues The 
new operators Pkyk^ automatically obey the fermionic 
anticommutation relationship Pk^k' } 

We get the ground state energy, 

2N 

{Hsp) = T.^k'^M-Eiykj+Y.^i 

7V3A2 



2[/, 



■(i + <5-q,q) 



(11) 



sp 



and the gap equation. 
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EEKt-^ikMEiykj- (12) 



Here, I labels the eigenenergy, and mi, m[ labels the 
matrix elements corresponding to the original S, P op- 
erators in the gap equation. 

The parameters used in the mean-field calculations are 
the same as in the ID case with small t-^^s added, and 
we still expect that the order parameter has the mo- 
mentum around 0.4357r/a as before. By self-consistently 
solving for Q and A, in the case = 0.05, the ground 
state is the cosine wave phase with Q = 0.4337r/a and 
A = 0.822. The ground state energy per site is —2.5927, 
lower than the noninteracting value —2.5896. When 
= 0.1, the ground state is also the cosine wave phase 
with Q = 0.4337r/a and A = 0.542. The ground state 
energy per site is —2.5955, lower than the noninteracting 
value —2.5949. These results confirm that (i) the cosine 
wave state has lower energy than the exponential wave 
state, (ii) the order parameter has the momentum close to 
the prediction of Eq. ([5]) , and (iii) larger transverse hop- 
ping tends to destroy the p-orbital pair condensate since 
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FIG. 2. (Color online) The occupation of s and p band 
within the paired state for different transverse hopping t±. 
Only the first quadrant of the Brillouin zone in the kx — ky 
plane is shown, kz = 7r/a. The black dashed lines indicate 
the "bare" Fermi surfaces for corresponding noninteracting 
fermions {Usp = 0). (a) (S+Sk) for = 0.05; (b) {SlSy,) for 
= 0.1; (c) (Pi^Pk> for = 0.05; (d) (Pj^Pk) for = 0.1. 



the energy gain for larger transverse hopping is mueh 
smaller than for smaller transverse hopping. 

An interesting feature of the p-orbital pair conden- 
sate in quasi-lD is the possible existence of Fermi sur- 
faces with gapless energy spectrum. We monitor the 
fermion occupation number, i.e. {Sl^S]t} and (-Pk^k) for 
increasing transverse hopping. The results are shown in 
Fig. m For small t^, they take the usual BCS form and 
vary smoothly from 1 (red) to (blue) across the bare 
Fermi surface (with interaction turned off), as shown in 
Figs. Ilja) and[5](c) for t-^ = 0.05. For larger transverse 
hopping, sharp Fermi surfaces characterized by a sud- 
den jump in (S'j[.S'k) and (^'^-Pk) appear. This is clearly 
shown in Figs. HKb) and Hfd) for t-^ = 0.1 as the oc- 
cupation number changes discontinuously from 1 (red) 
to (blue). It can be understood qualitatively as fol- 
lows. As i"*" increases, the original Fermi surfaces acquire 
a larger curvature in the transverse directions and the 
pairing condition in Eq. ([5]) cannot be satisfied every- 
where anymore. Therefore in some regions fermions are 
not paired and Fermi surfaces survive. One should also 
note that the calculation is based on the assumption that 
t-^ <^ t\ which predicts that Q is in the parallel direc- 
tion. This prediction should fail as t-^ increases beyond 
certain critical values. 



V. PHASE DIAGRAM 

Now, we systematically explore the phases of our sys- 
tem for general band filling and spin imbalance. Since 
we have s- and p- bands with different band widths, we 
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FIG. 3. (Color online) Band occupation for the four possi- 
ble phases in the system. The band colored in red represents 
the s band occupied by spin 4- fermions and the band colored 
in green represents the p band occupied by t fermions. The 
spin t fermions in the s band are not shown since they are 
inert, (a) Normal phase I (Nl) with one band empty and the 
other partially filled, (b) Normal phase II (N2) with one band 
fully filled and the other partially filled, (c) Commensurate 
p-orbital pair condensate (CpPC) with both bands partially 
filled. The occupation numbers are the same, (d) Incom- 
mensurate p-orbital pair condensate (IpPC) with both bands 
partially filled. The occupation numbers are different. 



introduce two dimensionless quantities for the chemical 
potentials fig and fip 

fis _ Ms 

2ts~ 2 ' 

fip ~ UJb _ — LUb 



2tp 16 



(13) 



Thus, for a non-interacting system, — 1 < < 1 con- 

trol the filling for the s and p-band fermions respectively. 
We then define the quantities 



h 
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(14) 



as the parameters controlling the average filling and po- 
larization in the phase diagram. The phase at — /i, —h is 
the same as the state at /i, h, since the transformation 
fi,h ^ —p,, —h gives fig, fip —fJ-s, — Mp, and the mean- 
field Hamiltonian with , fip is identical to Hamiltonian 
with —fis, —fJ-p via a particle-hole transformation up to a 
constant. 

We have four possible phases in such a system as shown 
in Fig. 131 As before, we ignored the inert fully filled s 
band of spin t fermions. We consider the p band of spin 
'I fermions and s band of spin ], fermions. When one of 
these two bands is empty and the other is filled, the pair- 
ing does not happen and we call it normal phase I (Nl) as 
in Fig. ^a). When one of these two bands is fully filled 
and the other is partially filled, the pairing also docs not 
happen since the fully filled band is inert. We call it nor- 
mal phase II (N2) as in Fig. |3jb). When both of them 
are partially filled, fermions near Fermi surfaces from the 
two bands will be paired and the system is in superfiuid 
phases as shown in Figs.[3l[c) and[3ljd). In the superfiuid 
regime, when h is small, the pairing momentum prefers 
Q = 7r/a and we call it commensurate p-orbital pair con- 
densate (CpPC). It is a special case of the p-orbital pair 
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condensate, where the occupation numbers of s-band spin 
], fermions and p-band spin t fermions are the same. It 
is similar to the conventional unpolarized pairing (BCS), 
where the spin '|~ fermions and spin ], fermions have the 
same population. However, in BCS pairing the CMM of 
the pair has the property Q = 0, while here Q ~ ixja. 
To understand the momentum tt/o preference, note that 
in conventional BCS case, the two species of fermions 
have the same energy spectrum and the pairing is be- 
tween two fermions with opposite momenta, which leads 
to the CMM of pair (3 = 0. Here, the structure of energy 
spectrum of p band is different from s band. The equal 
occupation numbers mean fci?-|-=7r/a — fc^j., which gives 
rise to Q = kp^ + kpi = ir/a, as shown in Fig. ^c). At 
last, when h is large, the pairing momentum stays at a 
general Q ~ kp-^+kp^ and the occupation number for the 
two species of fermions differ. We call it incommensurate 
p-orbital pair condensate (IpPC) as shown in Fig. ^d). 

To determine the phases, we minimize the free energy 
as a function of the pairing amplitude A and pairing mo- 
mentum Q by mean-field analysis using the cosine wave 
function as outlined in the previous section. When the 
minimum is realized at A = 0, it is normal phase. When 
A is finite, there are two possibilities. When Q = ir/a, it 
is CpPC. When Q n/a, it is IpPC. For the transition 
between superfluid and normal phase, and the transition 
between CpPC and IpPC, the behaviors of free energy 
show that the phase transitions are first order in a lat- 
tice system. Between the superfluid and normal phases, 
near the phase transition, A changes suddenly from 
to finite, and the free energy shows two local minima 
at A = and A 0. Between CpPC and IpPC, the 
pairing momentum changes from Q = n/a to Q ^ n/a 
discontinuously, and the free energy as a function of Q 
also has two local minima at Q = 7r/a and Q ^ 7r/a. 
Thus, they are first-order phase transitions according to 
our mean field analysis. Therefore, we can determine the 
phase boundaries between normal phase and superfiuid 
phase by monitoring A changing from zero to finite. We 
can also monitor Q changing from Q = Tr/a to Q^n/a 
to determine the phase boundaries between CpPC and 
IpPC. 

In Fig. m we present a phase diagram for = 0.05. 
The x's in Fig. 2] show the data points for the phase 
boundary obtained from the numerical procedure, and 
by connecting them we get the phase boundaries. An 
illustrative physical understanding about this phase di- 
agram is as follows. In Fig. SI when chemical potential 
difference h is small and the two bands are still partially 
filled to ensure the pairing, the system tends to stay in 
CpPC where Q = tx ja. It is similar to the conventional 
BCS superfiuid case. As h becomes larger, as long as the 
average filling \i is not too large or small and the two 
bands are still both partially filled, the pairing persists 
despite the spin imbalance and the system is in IpPC. 
If \i gets more and more negative, the average filling be- 
comes smaller and smaller, and at certain /i, /i, -p band 
of spin t fermions will be empty and the system will be- 
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FIG. 4. (Color online) The phase diagram of the p-orbital pair 
condensate for = 0.05. /i and h are defined in the main 
text. CpPC: the s band of spin \. fermions and the p band of 
spin t fermions have the same occupation numbers. IpPC: the 
s band of spin \. fermions and the p band of spin f fermions 
have different occupation numbers. Nl with the p band of spin 
t fermions empty and the s band of spin 4- fermions partially 
filled. N2 with the p band of spin t fermions partially filled 
and the s band of spin \. fermions fully filled. 

come Nl without pairing. Similarly, when [i is large and 
positive, the average filling is very high and at certain 
/i, /i, the s band of spin 4- fermions will be fully occupied, 
and the system becomes N2 without pairing. The almost 
straight phase boundaries in Fig. U between IpPC and 
normal phases indicate that these phase transitions are 
due to the change of band occupation as empty o par- 
tially filled o fully filled. In Fig. |4l the phase boundary 
between IpPC and Nl corresponds to the critical condi- 
tion that the s band of spin \, fermions is partially filled 
while the p band of spin f fermion becomes empty, and 
the almost straight phase boundary corresponds to the 
condition that p,p — j.i — h ^ —I (but, as before, this 
is only an approximate argument due to the presence of 
interaction) . Similarly, the almost straight phase bound- 
ary between IpPC and N2 corresponds to the condition 
that the s band of spin J, fermions becomes fully filled, 
while the p band of spin f fermions is partially filled, or 
fig = fJ.+h = 1. All the phase transition lines in Fig.|3]are 
mean field results, and these straight lines are expected 
to be corrected by quantum critical fiuctuations. The 
phase diagram shows that the p-orbital pair condensate 
happens in large parameter regimes and is closely related 
to the band and orbital properties in the optical lattice 
systems. 



VI. SIGNATURE OF THE p-ORBITAL PAIR 
CONDENSATE IN MOLECULE PROJECTION 
EXPERIMENT 

The p-orbital pair condensate phase can inspire impor- 
tant experimental signatures for finite momentum con- 
densation of bosonic molecules in higher orbital bands. 
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FIG. 5. (Color online) (a) The momentum distribution func- 
tion Uq of projected molecules for a quasi-lD system with 
= 0.05 (all other parameters are same as before) accord- 
ing to mean field theory. Here, q — q^, Qy — Qz = 0. 
(b) Pair correlation function Cq for a ID chain of A'^ = 60 
sites obtained by DMRG. The peak is located at 0.4337r/a 
in both figures, which corresponds to the value kpi^ + kpi = 
[Ns+N - Np)n/Na for A^, = 49 and Np = 15. 



By fast sweeping the magnetic field (and thus the interac- 
tion) from the BCS region to the deep EEC region across 
a Feshbach resonance, the ECS pairs arc projected onto 
Feshbach molecules, which can be further probed for ex- 
ample by time-of-flight images [13| . The bosons produced 
effectively reside in p band and are stable, since by Pauli 



blocking the filled s-band fermions will prevent the the 
p-wave bosons from decaying [l^ . Here, we use a simple 
model [2^ to evaluate the momentum distribution of 
molecules after projection 

= X! /k/k'('5'k+q/2-P-k+q/2-P-k'+q/2S'k'+q/2), (15) 
k,k' 

where /k is the molecular wave function, and the correla- 
tion function can be evaluated within mean field theory 
[26| . For fast sweeps, the molecular size is small com- 
pared to lattice constant and its wave function can be 
approximated by a delta function in real space (a con- 
stant \Jl/N in momentum space). Ey this assumption, 
Uq is the same quantity as Cq in Eq. (|6]). Figure [5l^a) 
shows the rig of p-wave Feshbach molecules and a peak is 
located at 0.4337r/a. FigurejSJb) shows Cq from Eq. 
based on the DMRG results shown in Fig.[lJb).The time- 
of-flight experiment is predicted to distribute peaks cor- 
responding to that in Fig. [51 Note that for the ID prob- 
lem (Fig. [5{b)), the delta- function peak is replaced by a 
cusp characteristic of power law due to the lack of long 
range order. 

We thank Chungwei Lin for helpful discussions. This 
work is supported by ARO Grant No. W911NF-07-1- 
0293. 
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